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ABSTRACT 

We investigate the presence of iron line emission among faint X-ray sources identified in the IMs Chandra Deep 
Field South and in the 2Ms Chandra Deep Field North. Individual source spectra are stacked in seven redshift 
bins over the range z = 0.5-4. We find that iron line emission is an ubiquitous property of X-ray sources up to 
z~3. The measured line strengths are in good agreement with those expected by simple pie-Chandra estimates 
based on X-ray background synthesis models. The average rest frame equivalent width of the iron line does not 
show significant changes with redshift. 

Subject headings: Surveys - Galaxies: active — X-rays: galaxies -X-rays: general -X-rays:diffuse background 



1. INTRODUCTION 

It has been pointed out already a decade ago (Matt & Fabian 
1994) that prominent spectral features, and especially the ~ 6.4 
keV iron Ka emission line commonly found in Seyferts spec- 
tra (Nandra & Pounds 1994), may lead AGN synthesis models 
to predict a detectable signature in the spectrum of the cosmic 
X-ray background (XRB) around a few keV. 
The integrated contribution of iron line emission in the sources 
making the XRB has been quantitatively estimated by Gilli et 
al. (1999), on the basis of a large database of observational re- 
sults available in the pre-Chandra and XMM-Newton era. The 
predicted shape of the integrated emission from individual lines 
was computed following the prescriptions of the XRB synthe- 
sis model of Comastri et al. (1995) resulting in a rather broad 
bump extending from about 1 keV (depending on the maximum 
redshift assumed) up to 6.4 keV. The expected intensity of such 
a bump above the XRB level is below the 5% level over most 
of the energy range and reaches a maximum value of the or- 
der of 7% around 2-3 keV, the exact value being dependent on 
the redshift (Zcm) at which the evolution of the luminosity func- 
tion is stopped. Although such an estimate is model dependent, 
Gilli et al. (1999) concluded that it can be safely regarded as an 
upper limit to the total iron line contribution. 
In the last few years our knowledge of the XRB sources has 
been significantly improved by deep Chandra and XMM- 
Newton surveys. A fraction as large as 80-90% of the XRB 
flux below 5-6 keV has been resolved into individual sources 
(Worsley et al. 2004). Extensive multiwavelength follow-up 
observations have clearly established that the X-ray source red- 
shift distribution is peaked at z ~0.7-l (e.g. Hasinger 2004, 
GiUi 2004) and that about 60% of the XRB originates at z < 1 
(Barger et al. 2003). In the XRB synthesis model worked out 



by Gandhi &. Fabian (2003) the iron line contribution turned 
out to be maximum at about 3.8 keV (corresponding to a typi- 
cal redshift of 0.7) and the excess due to the iron line is of the 
order of 3^ %. 

Despite the increasing number of independent observations 
of the XRB spectrum below 10 keV with BeppoSAX (Vecchi 
et al. 1999), XMM-Newton (De Luca & Molendi 2004) and 
RXTE (Revnivtsev et al. 2004), yielding good signal to noise 
ratio measurements of the extragalactic background, the accu- 
racy reached so far is not such to detect the iron features at the 
level expected by the model predictions outlined above. 

An alternative approach, devised to avoid the line smearing 
due to the large redshift range over which AGN spectra are 
summed and the present uncertainties in the XRB spectrum, 
is to search for iron features over appropriate redshift bins. The 
detection of an iron line and the study of its intensity and pro- 
file would open the possibility to investigate interesting issues 
such as the metal abundances and the relative fraction of Comp- 
ton thick sources (which are known to have extremely strong 
iron lines) at high redshift, along with relativistically broadened 
lines beyond the local Universe (Comastri, Brusa, Civano 2004; 
Streblyanska et al. 2005). In the following we will present and 
discuss the results obtained for a large sample of spectroscop- 
ically identified sources in the X-ray surveys which resolved 
most of the XRB: the 2 Ms pointing in the Chandra Deep Field 
North (Alexander et al. 2003; hereinafter CDFN) and the 1 Ms 
observation of the Chandra Deep Field South (Giacconi et al. 
2002; hereinafter CDFS). Throughout the paper, the adopted 
values for the Hubble constant and the cosmological parame- 
ters are Ho=70 km s"' Mpc"', rjA=0.7, O,n=0.3. 

2. DATA ANALYSIS 
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2.1. The sample 

The optical spectroscopic and photometric redshifts of the 
X-ray sources detected in the CDFN are presented by Barger et 
al. (2003), while a list of spectroscopic redshifts of the sources 
detected in the CDFS is reported in Szokoly et al. (2004; 
hereinafter S04). Photometric redshifts based on good qual- 
ity HST/ACS images and deep ISAACA^LT observations are 
available for all but four CDFS sources without spectroscopic 
redshift (Zheng et al. 2004; hereinafter Z04), leading to a vir- 
tually complete catalog of identified X-ray sources. Here we 
have considered only those CDFS photo-z with a quality flag 
greater than 0.5 (see Z04 for details) ensuring a reliable red- 
shift estimate. Moreover, given that both S04 and Z04 identi- 
fications are based on the X-ray source catalogue of Giacconi 
et al. (2002), we have revised some optical/X-ray associations 
according to the improved astro metry provided by Alexander et 
al. (2003). Only sources with spectroscopic redshifts have been 
included in the CDFN sample because no details are given on 
the photo-z quality in Barger et al. (2003). The larger CDFS 
redshift completeness allows us to obtain a counting statistic 
per redshift bin similar to that of the CDFN despite the shorter 
exposure. 

Individual source spectra have been stacked together in seven 
redshift bins spanning the z=0.5^ range (Table 1). The choice 
of bin sizes and distribution is driven by a trade-off between 
the number of counts in each bin and the need to sample a cor- 
responding observed energy range narrow enough to detect the 
spectral feature, keeping at the same time the instrumental re- 
sponse as uniform as possible. For these reasons, the redshift 
interval below z=0.5 is not considered; moreover, the redshift 
interval 1.7 < z < 2.5 is excluded because the iron line is red- 
shifted in the ^ 1 .8-2.4 keV energy range, which encompasses 
the sharp drop in the effective area due to the instrumental irid- 
ium edge. 

The iron line signal is not expected to be significantly diluted 
by non-AGN sources. Below z=2, the Ka Une faUs in the hard 
band (> 2 keV), where the detected sources are mostly AGN. 
At higher redshifts it is shifted to the soft band, but then only 
sources with AGN luminosities can be detected. In order to 
avoid contamination by a few individual bright sources which 
could dominate the stacked signal we have excluded from each 
bin the brightest object, if its contribution exceeds 50% of the 
entire flux in the bin. The chosen redshift bins, the expected 
position of the 6.4 keV line, the bin width and the number of 
sources in each bin are listed in Table 1 . The final sample in- 
cludes 171 sources in the CDFN and 181 in the CDFS, spanning 
the luminosity range Li-iofev = 10"*' - 10'*^ erg s"' (see Table 1). 
The small solid angle (~ 0.2 deg"^) covered by the deep fields 
does not allow us to sample the AGN population at higher lu- 
minosities. 

' http://cxc.harvard.edu/ciao/ 

^ http://cxc.harvard.edU/ciao3.0/threads/merge_all/ 

' http://hea-www.harvard.edu/alexey/acis/tgain 

We have also verified that applying the new gain correction tool released with the CALDB 2.28 version the results are not modified. 
' See also: http://cxc.harvard.edu/ccw/proceedings/03_proc/presentations/marshaI12 



2.2. X-ray stacking 

The X-ray data for both the CDFN and CDFS observations 
have been retrieved from the public archive and processed with 
standard tools making use of the calibrations associated with 

the CIAO' software (version 3.0). 

A total of 20 and 1 1 pointings in the CDFN and CDFS, respec- 
tively, were summed together with the merge_all^ script in 
two different merged event files. The Charge Transfer Ineffi- 
ciency (CTl) and gain^ corrections were also applied to each 
single pointing"^. Spectra, response matrices and effective ar- 
eas of the stacked spectrum in each redshift bin were extracted 
from the merged event files using the standard CIAO tools de- 
veloped to properly weight responses and effective area files for 
multiple extraction regions (e.g. mkwar f ; see Civano, Comas- 
tri & Brusa 2005 for a detailed description of the procedure). 
The target positions in each redshift bin span almost all the of- 
faxis angles, thus the extraction radius was varied from 4 to 16 
arcsec depending on the source brightness and off-axis posi- 
tion. Several background spectra were extracted with the same 
procedure from a stack of regions nearby each source in each 
redshift bin and varied in size and shape. 
Finally, the stacked spectra in each redshift bin of both 
the CDFN and CDFS were summed together with stan- 
dard FTOOLS routines (mathpha, addarf, addrmf) 
weighted for both the exposure time and the number of 
counts. In the foUowing we will consider only the resulting 
CDFN-hCDFS spectra. 

3. SPECTRAL ANALYSIS 

Since we are interested in emission features over a well de- 
fined energy range, source spectra have been extracted in the 
1-6 keV band, to minimize calibration uncertainties in the low 
energy response. X-ray spectra were rebinned to have at least 
20 counts per bin and fitted with XSPEC (version 11.3.0; Ar- 
naud 1996); errors are reported at the 90% confidence level for 
one interesting parameter (Ax^=2.71). 

The continuum is parameterized with a single power-law with 
the slope free to vary. A negative (r = -0.17) edge at 2.07 keV 
has been added to account for residual uncertainties in the effec- 
tive area above 2 keV, as suggested by the CXC team (VikhUnin 
et al. 2005)^. With the exception of the highest redshift bin, a 
significant excess above a power-law continuum is present over 
the energy range 6A/(l+Znmx) - 6A/(l+Zmin) keV, where Zmin 
and Zmax are the bin boundaries (see Fig. 1). We then added a 
redshifted gaussian component (model po+zgauss) with the 
width free to vary, centroid energy fixed at 6.4 keV and redshift 
fixed at the mean value of each bin (z), and verified that in all 
cases the fit improved, though with a different level of signifi- 
cance. Table 1 lists the best fit parameters for each redshift bin. 
The equivalent widths (EWs) quoted in Table 1 are measured 
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in the observed frame. To get the intrinsic, rest frame EWs a 
multiplicative factor 1 -I- z must be applied*. 
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3.1. Safety checks 

In order to check whether the results are dependent from the 
counting statistics and/or the energy range over which the con- 
tinuum underlying the iron excess is measured, we have per- 
formed several safety checks. 

First, we have further excluded in each bin the two or three 
brightest sources, if their contribution to the total counts was 
larger than 60%. Then, spectral fits were repeated in the 2.2- 
6 keV energy range for sources at z < 2. Finally, a redshifted 
absorption component (assumed to be at z in each redshift bin) 
was added to the continuum. In all cases, the spectral parame- 
ters were found in agreement, within the statistical errors, with 
those reported in Table 1. Furthermore, for the z< 2 bins for 
which the statistics is higher, we verified that, when the redshift 
of the line component is left free to vary, the resulting best fit 
redshift is always almost coincident with that expected from a 
line at 6.4 keV. 

Extensive simulations have also been performed to quantita- 
tively assess the line broadening introduced by redshift smear- 
ing. The spectra of the sources contributing to each bin were 
modeled as a single power law plus a narrow unresolved (cr=0. 1 
keV) Ka line. The actually observed sources redshift distribu- 
tion within the bin is further assumed and the individual sim- 
ulated spectra are stacked together. With the exception of the 
z=0. 9-1.1 bin, all the measured widths are in agreement with 
those resulting from the simulations, and therefore consistent 
with being produced by the superposition of narrow features 
(but see also the discussion). 

According to XRB synthesis models, a large fraction of the 
sources contributing to each redshift bin are obscured by col- 
umn densities in the range r^lQ^^~^'^ cm"^. When obscured 
sources at different redshifts are stacked together, the iron line 
underlying continuum is expected to be modified, with respect 
to a single power law, by the most prominent absorption fea- 
tures such as the low energy cut-off and the 7.1 keV iron edge. 
To take this into account, we assumed a continuum shape as re- 
sulting from the Gilli et al. (2001) model once the AGN lumi- 
nosity function is integrated in the narrow redshift bins adopted 
here. Although the recomputed EW are lower than those re- 
ported in Table 1 (by about 30-50%), the line excess is still 
significant. A more detailed investigation would require an ex- 
tensive analysis of the AGN synthesis models parameter space 
and is beyond the purposes of this letter. 

* The superposition of sources at different redsliift in eacli bin should return a slightly 
by using (1 + z) is negligible. 
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Fig. 1. — Residuals of a simple power-law fit to the source spectra in 
seven different redshift bins as labeled. The vertical line in each panel is at 
the expected position for the redshifted 6.4 keV Fe Ka line while the shaded 
region encompasses the bin width reported in Table 1 and defined as AE = 
6.4/(l+z„,„,-)-6.4/(l+z„„„)keV. 



4. DISCUSSION 

The rest frame equivalent widths for the seven redshift bins 
are compared (Fig. 2) with those predicted by the XRB syn- 
thesis model of Gilli et al. (2001, shaded area in Fig. 2). The 
line EW distribution is modeled following Gilli et al. (1999, 
see their Table 4) and is a function of the absorption column 
density increasing from about 250 eV for unobscured sources, 
to 400 eV for obscured (logA^n < 24) AGN, reaching about 2 
keV for Compton thick sources. These values are constant over 
the entire redshift range. The upper bound of the shaded re- 
gion has been computed assuming EW=280 eV for unobscured 
sources, consistent with the average value measured by ASCA 
in the nearby Universe. Recent XMM-Newton and Chandra 
results have questioned the presence of broad iron lines in 
many unobscured-type 1 AGN, measuring in turn a lower EW 
(^ 100- 150 eV on average) along with a decrease of the line 
EW towards high luminosities (e.g. Page et al. 2004). We 
have quantitatively estimated the effects of reducing the line 
EW among unobscured sources by considering the most con- 
servative case of no iron emission at all (lower bound in Fig. 2). 
Since in our model most of the iron signal is produced by ob- 
scured sources with 22 < log(NrmH)< 24 cm"^, the overall EW 
are reduced by only 15 - 20%. 

different multiplicative factor, but the difference in the rest frame EW obtained 
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Fig. 2. — Rest frame equivalent widths as a function of redshift compared 
witli the model predictions (shaded strip; see text for details). 

The lack of any dependence of the average line intensity 
upon redshift up to z ^ 3 could be naively interpreted as a con- 
stant iron abundance with redshift. The real case is likely to 
be much more complicated, given the several factors affecting 
the detected signal. As an example, the population of heavily 
obscured Compton thick sources (which have the strongest iron 
lines) might be hiding among the X-ray sources either unidenti- 
fied or with a poorly constrained photo-z and thus not included 
in the present analysis. However, according to our calculations 
(Gilli et al. 2001), Compton Thick AGN are expected to pro- 
vide only about 5-10% of the iron signal in each redshift bin. 
The presence of a population of Compton Thick sources much 
larger than that assumed would be more easily appreciated by 
the distortions introduced on the broad band 3-100 keV XRB 
spectrum rather than by a stronger line in the various redshift 
bins. 

A significant decrease of the line EW towards high luminosi- 
ties has been reported by Nandra et al. (1997) for a sample 
of type 1 AGN observed with ASCA. The effect is more pro- 
nunced for luminosities higher than 10"*^ erg s"' . Unfortunately, 
our sample is not well suited to investigate iron line proper- 
ties at high luminosities due to the lack of sources luminous 
enough (see Table 1) to check the Nandra et al. (1997) find- 
ings. Moreover, the stacked spectra are dominated by obscured 
sources rather than type 1 AGN adding further uncertainties in 
the comparison. 



Our findings are also consistent with those reported by Stre- 
blyanska et al. (2005) from the analysis of the rest-frame 
stacked spectrum of identified sources in the XMM-A^ewfo« ob- 
servation of the Lockman Hole (~ 400 eV, when a gaussian pro- 
file is considered). The residuals leftward of the iron line which 
in our fits are more prominent in the bins with the highest num- 
ber of sources and counting statistic (z=0.5-0.7 and 0.9-1.1) 
suggest the presence of a broad redshifted component, similar 
to that observed by Streblyanska et al. (2005). Since the main 
goal of the present analysis is the detection of iron emission in 
the sources of the XRB up to high redshifts, our approach is 
not designed to investigate in detail the average line shape. We 
recall, however, that the observed EW and line profile are de- 
pendent on an accurate modeling of the underlying continuum. 
The superposition of absorbed spectra with different redshifts 
has the effect of reducing the measured EW (see previous Sec- 
tion) and furthermore may produce a spurious red wing. 

5. SUMMARY 

We report the detection of iron emission up to z ^ 3 in the 
X-ray spectra of faint Chandra sources stacked into different 
redshift bins. The measured EW are in agreement with a sce- 
nario in which the lines are intrinsically narrow and their inten- 
sity does not change significantly with redshift (and/or luminos- 
ity), which can be interpreted as a constant iron abundance as 
a function of redshift. Extensive simulations and safety checks 
have been performed in order to test the reliability of our re- 
sults. Even taking into account the effects of absorption fea- 
tures which modify the high energy power law in obscured 
sources the detection of the emission line remains significant. 
Despite their prominent iron lines, our approach does not al- 
low to put tight constraints on the number density of Compton 
Thick sources, that suffer from extreme absorption in the Chan- 
dra band and can be detected only in small numbers compared 
to Compton thin AGN. Although there might be hints for the 
presence of gravitationally redshifted broad line components, 
we caution that their intensity and profile significantly depends 
on the modeling of the underlying continuum. 
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TABLE 1 

Redshift bins and number of sources 



min-max 


E(Kaf^> 
keV 


Bin width'-'' 
keV 


Number of sources'^* 
CDFS CDFN Total 


lO**^ erg,s-' 


Counts'"' 




PL -1- zgauss'"' 
X^/d.o.f. 


Ax^'« 


^(10) 
keV 


EW*"> 
eV 


0.5-0.7 
0.7-0.9 
0.9-1.1 
1.1-1.4 
1.4-1.7 


4.00 
3.55 
3.20 
2.78 
2.56 


0.50 
0.40 
0.32 
0.38 
0.30 


45(7) 
34(10) 
30(11) 
25(10) 
19(11) 


39 
32 
58 
25 
4 


84 
66 
88 
50 
23 


1(0.1-140) 
4(0.3-300) 
8(0.6-400) 
20(0.9-400) 
40(2-500) 


16831 
13240 
34958 
15629 
7841 


1.20±0.03 
1.44±0.04 
1.47±0.02 
1.45±0.04 
1.36±0.04 


284.2/289 
229.7/260 
342.6/321 
258.3/267 
248.0/195 


44.2 
14.0 
112.6 
63.8 
19.0 


0.27!S;S 
0.48*1^ 
0.3ltl^ 
0.18t|^ 


207!55 

116_58 

223*^^ 

101!^^ 


2.5-3.0 
3.0-4.0 


1.71 
1.42 


0.23 
0.32 


18(9) 
9(3) 


4 
9 


22 
18 


100(20-1000) 
200(36-2000) 


4895 
2919 


1.46±0.07 
1.30±0.10 


164.5/141 
82.3/1 1 1 


10.0 
2.0 


'■'-'-0,40 

< 1.94 


< 89 



Notes: '''Redshift range for the stacked spectra. *^'Expected position of the Fe 6.4 keV line in the observed frame. *^*Observed frame bin width, in keV, at the 6.4 keV position. '"^'Number 
of sources in the CDFS, CDFN and combined samples (in parenthesis for the CDFS sources with photometric redshifts). ^^^Median 2-10 keV luminosity (in parenthesis the minimimi and 

maximum luminosities) '^^Total net counts in the 1-6 keV band. '^'Best fit power-law spectral index and 90% errors. '^'Statistical significance of the spectral fit for the po+zgauss model. 
Ax^ with respect to the single power-law model. ^^'^'^^'Best fit line width (a) and observed frame intensity (EW) of the gaussian component and 90% errors. 



